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ABSTRACT 



We present a detailed analysis of a long, 100 ksec, observation of NGC 4151 by ASCA, con- 
temporaneous with an observation by the CGRO/OSSE. We fit the data with physical models 
including an Fe K line and Compton reflection both relativistically broadened and coupled ac- 
cording to theoretical results. The model also includes a narrow Fe K component, which is 
emitted by both an extended plasma region and the X-ray absorber. Our study of the absorber 
shows strong evidence for the presence of more than one partial-covering cloud in the line of sight. 
Taking these points into account, our best intrinsic model includes a Comptonization continuum 
with an X-ray slope of F ~ 1.9 from a thermal plasma with an electron kT ~ 70 keV and a 
disk line with an equivalent width of ~ 70 eV. The broadening of the line and reflection indicate 
their origin from innermost parts of an accretion disk. Our results indicate that NGC 4151 has 
an intrinsic X-ray spectrum and variability properties similar to those of average Seyfert Is. 
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NGC 4151, a Seyfert 1 galaxy at z = 0.0033, 
has the ~ 2-10 keV photon index of its power- 
law component reported from observations by sev- 
eral different satellites to be in the range of 0.3 < 
r < 1.8 (Yaqoob & Warwick 1991; Yaqoob et al. 
1993; Wang et al. 1999, hereafter W99; Ogle et al. 
2000; Yang et al. 2001; Piro et al. 2002). (Note 
that the lowest of these values, 0.3-0.4, have been 
measured by Chandra.) The corresponding aver- 
age Seyfert 1 photon index is F = 1.95 with an 
intrinsic standard deviation of only 0.15 (Nandra 
& Pounds 1994), so the lower values in the above 
range are quite unusual. Since NGC 4151 is the 



^Centrum Astronomiczne im. M. Kopernika, Bartycka 
18, 00-716 Warszawa, Poland; aaz@camk.edu.pl 

■^Department of Physics and Astronomy, The University 
of Oklahoma, 440 W. Brooks St., Norman, OK 73019, USA; 
leighly@ou.edu 

^NASDA-SURP 2-1-1, Sengen, Tsukuba, Ibaraki 305- 
8505, Japan 

■^Istituto Te.S.R.E, CNR, Via Gobetti 101, 40129 
Bologna, Italy 

5RIKEN, 2-1 Hirosawa Wako, Saitama, 351-0198, Japan 



most often studied Seyfert (as it is the brightest 
one in the hard X-rays), it is of great interest to 
determine whether it is an archetype of its class 
or very far from it, as would be indicated by the 
spectral index varying up to llcr from the average. 

Reliable determination of the spectral indices 
of power law indices in accreting black holes is 
of major theoretical importance. Our present un- 
derstanding of the power law component in both 
Seyferts and black hole binaries associates its ori- 
gin with the process of thermal Comptonization 
(e.g., Poutanen 1998; Zdziarski 2000). However, 
this process yields in general F > 1 (e.g., Sun- 
yaev & Titarchuk 1980). If, on the other hand, 
Comptonization were non-thermal, F > 1.5 (in- 
cluding cooling and pair production, e.g.. Light- 
man & Zdziarski 1987). Very hard local continua 
can be obtained from Comptonization in a highly 
relativistic plasma {kT > 500 keV) of a small op- 
tical depth (e.g., Haardt 1993), but this model for 
NGC 4151 is then in conflict with the measure- 
ment of its plasma temperature of kT < 100 keV 
(Johnson et al. 1997, hereafter J97; Zdziarski et 
al. 2000). If the measurement of F ~ 0.3-0.4 in 
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NGC4151 (Ogle ct al. 2000; Yang ct al. 2001) is 
indeed confirmed, it would call for a new physi- 
cal explanation. Furthermore, Yang et al. (2001) 
found no contribution of the scattered component 
to the soft excess in NGC4151. However, this re- 
sult was based on their value of F = 0.3 obtained 
assuming a single uniform absorber. On the other 
hand, the fitted value of F strongly depends on the 
model of the complex absorber in this source, as 
noted by Weaver et al. (1994, hereafter W94). 

Also, the Fe Ka line has been found to be 
very broad and strong in ASCA observations 
during hard states of NGC4151 (Yaqoob et al. 
1995; W99). A small F during those observations 
would constitute an interesting counterexample to 
the positive correlation between the strength and 
relativistic broadening of Fe K reprocessing fea- 
tures and F observed in Seyferts (Zdziarski et al. 
1999; Lubinski & Zdziarski 2001), which is both 
expected theoretically and confirmed to hold in 
black- hole binaries (Gilfanov et al. 2000). 

On the other hand, monitoring by the OSSE 
detector aboard CGRO shows the shape of its 
50-300 keV spectrum to be weakly variable and 
very similar to the corresponding average spectra 
of Seyferts (J97; Zdziarski et al. 2000). In par- 
ticular, the model of thermal Comptonization and 
Compton reflection (standard Seyfert continuum 
components) applied to the average OSSE spec- 
trum of NGC4151 yields the 2-10 keV index of 
F ~ 1.9 (J97), i.e., very close to the measured av- 
erage of Seyfert Is. 

It is quite possible that the unusual proper- 
ties of NGC4151 noted above arise because the 
X-ray spectrum of NGC4151 is usually strongly 
absorbed with a characteristic column density of 
A^H ~ 10^^ cm~^, which makes measuring the in- 
trinsic X-ray index difficult (as already noted by 
W94). For instance, the observation with the his- 
torical weakest absorption (Yaqoob et al. 1993), 
A^H ^ 10^^ cm~^, has also the highest historical 
2-10 keV flux, a soft F ~ 1.8 and a low upper limit 
on the Fe K line (see also Zdziarski et al. 1996, 
hereafter Z96). This may indicate that an appar- 
ent correlation between the 2-10 keV flux and F 
(e.g., Yaqoob et al. 1993) is to a large degree due 
to strong absorption, which both suppresses the 
emitted X-ray flux and hardens the spectrum (al- 
though Yaqoob & Warwick 1991 have rejected that 
possibility based on their treatment of the absorp- 



tion). 

Z96 have suggested previously that the low val- 
ues of F and/or the very strong Fe K line may be 
artifacts of insufficient account of complex absorp- 
tion and neglect of Compton reflection. However, 
as noted by W99, the data used by Z96 were either 
of a low signal-to-noise ratio or limited energy res- 
olution. Subsequently, W99 found the Fe K line to 
be very strong and broad and F ~ 1.3 in a long, 
100 ksec, ASCA observation in 1995, and com- 
mented on the unprecedented high signal to noise 
of that observation. 

Here, we re-examine the above issues by rean- 
alyzing the 1995 ASCA observation. Our treat- 
ment is different from that of W99, who considered 
only the 1-10 keV SIS data, as we analyze the full 
ASCA spectrum together with a contemporaneous 
OSSE observation and employ relatively more so- 
phisticated physical models. Preliminary results 
of our study of the ASCA spectrum were given in 
Leighly et al. (1997, 1998). In our present best 
model, the intrinsic F ~ 1.9, equal to that ob- 
tained from modelling the soft 7-rays. We con- 
clude that previous reports of very low values of 
F in NGC4151 are likely to be due to simphfied 
treatments of the X-ray absorption and neglect of 
Gompton reflection. We find the broad line to have 
a moderate equivalent width of ~ 10^ eV, which 
is average for Seyfert Is with that F (Lubinski & 
Zdziarski 2001). Also, we find that the energy 
dependence of the variability amplitude implies 
a major role of complex variable absorption, in 
agreement with the results of the spectral fitting. 

2. Data Analysis 

NGG4151 was observed by ASCA from 1995 
May 10 02:49:43 tiU 1995 May 12 13:50:36 (UT, for 
the SIS) with the net exposure of ~ 10^ s for each 
of the detectors. SIS 1-CCD faint mode was used 
and standard data reduction was applied. Gare 
was taken to exclude photons from the BL Lac 
object 1E1207.9+3945 (which lies 5' from NGG 
4151), as shown in Figure 1. We have rebinned 
the original channels to 1 /4 of the FWHM resolu- 
tion of the SIS and GIS detectors, and assumed a 
1.5% systematic error, which was estimated from 
residual discrepancies among the 4 ASCA detec- 
tors. We note that this observation was made only 
2 years after the ASCA launch, and thus was only 
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Fig. 1.— The part of the ASCA GIS2 image of 
the NGC 4151 used in the present reduction. The 
region around the nearby BL Lac 1E1207. 9+3945 
has been removed. 

slightly affected by the time-dependent degrada- 
tion of the SIS resolution. 

We fit the data using XSPEC (Arnaud 1996). 

We assume elemental abundances from Anders & 
Ebihara (1982) and the interstellar Galactic ab- 
sorbing column of 2.17 x 10^° cm~^ (Murphy et 
al. 1996). The single-parameter uncertainties be- 
low correspond to 90% confidence (Ax^ = 2.7). 

3. The Average Spectrum 

3.1. Simple phenomenological models 

We start with a model similar to one of those 
used by W94, consisting of a power law, assumed 
here to be e-folded with = 400 keV, absorbed 
uniformly with a column density, Aq, and partially 
with a column A'l covering a fraction /i (the so- 
called dual absorber). The soft component is mod- 
eled by the sum of a fraction, /sc , of the power law 
scattering towards the observer (suffering only the 
Galactic absorption) and thermal bremsstrahlung. 
The model also includes a Gaussian Fe Ka line ab- 
sorbed only by the Galactic A^h, which is fitted to 



be narrow with the width of (t„ ~ 0.1 keV and the 
peak energy of ~ 6.3 keV. 

Such a model provided a good fit {x^ /y ~ 1) 
to the two ^ 20 ksec net-exposure data taken in 
1993 (W94). However, it provides a very poor 
model to the present data, /y = 1016/641, and 
it yields strong localized residuals, including one 
around ~ 5-6.5 keV, as illustrated in Figure 2a. A 
likely reason for this difference is the much better 
statistics of the present data, which has a longer 
exposure by a factor of > 5. 

In particular, including (in addition to the nar- 
row Gaussian) a broad Gaussian feature (absorbed 
as the power law) with the peak energy at i?b = 
5.9+[J;^ keV, the width of = 0.9toi keV, and 
equivalent width of Wb = 390 eV, reduces by 
141. Furthermore, the data show line- like resid- 
uals at 0.89 keV, and addition of a narrow line 
further reduces x^ by 123. This feature probably 
corresponds to the O vili radiative recombination 
continuum and Ne ix emission that originates in 
the X-ray emitting narrow-line region, and is seen 
clearly in the Chandra spectrum (Ogle ct al. 2000). 
Also, we find that the SIS continuum is system- 
atically harder than the GIS one. This is likely 
to be due to the residual inaccuracy of the SIS 
calibration for the present very hard X-ray spec- 
trum, and is recognized here because of the very 
good statistics. Addition of a small power law 
correction to the SIS model, AF, reduces x^ by 
41 (chance probability is 2 x 10~^; hereafter, we 
use the F-test to estimate significance of adding 
model components). Hereafter, we keep the soft 
line energy and AF fixed at 0.89 keV and -t-0.047, 
respectively. We note that W99; who analyzed 
the SIS data at energies > 1 keV only, have not 
considered the complexities of the soft spectrum 
present in the data. 

This model, with x^/^ = 711/637, which is 
within ~ 2(7 of the mean of the x^ distribu- 
tion, is formally acceptable but it has a num- 
ber of problems. First, at the bremsstrahlung 
temperature of 0.31 keV, strong line emission is 
expected. However, using the mekal model in 
XSPEC (Mewe et al. 1985) actually worsens the 
fit by Ax^ = +16. Second, the origin of the 
very strong and broad Gaussian line is unclear. 
It is indicative of very strong Compton reflection 
(e.g., Zycki & Czerny 1994), which has not been 
included in the model. The continuum spectral 
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Fig. 2. — (a) The residuals to the model with a power law, dual absorber, a Gaussian line and bremsstrahlung. 
(b) The residual for our best model for the ASCA data, obtained by adding a disk line with associated 
reflection, ionization of the full screen, and an additional partially-covering cloud. (The assumed systematic 
error of 1.5% is not shown.) 



index is F = LSO^q obi Compton reflection 
is negligible in other Scyfcrts with a such a hard 
spectrum (Zdziarski et al. 1999). Thus, this model 
can be considered only an empirical (rather than 
physical) description of the spectrum. 

3.2. Fe K line profile and Compton reflec- 
tion 

To make further progress, wo havc^ tric^l to de- 
velop a physical description of the spectrum by 
taking into account the following issues. The for- 
mation of Fc fluorescence lines proceeds through 
photoionization, so that the number of photons 
removed beyond the iron K-edge times the fluo- 
rescence yield equals the number of line photons. 
Broad Fe Ka lines are most likely formed by re- 
flection from a cold medium (Nandra & Pounds 
1994). In this case, basic physics demands that the 
line flux should be tied to the relative strength of 
Compton reflection, f2/27r (with the coefRcient of 
proportionality dependent on the ionization state; 
George & Fabian 1991; Zycki & Czerny 1994; 
Zycki et al. 1998). Furthermore, if a line is broad- 
ened and shifted kinematically and by gravity, the 
same broadening will occur to the reflection com- 
ponent, in particular its K edge. Then, mod- 
els with broad lines originating from illumination 



of an accretion disk but either without reflection 
or with static reflection are not physically self- 
consistent. 

We note that there is evidence that in some 
Seyferts with observed broad Fe K lines the sim- 
ple model with irradiation of an accretion disk ex- 
tending close to the minimum stable orbit may 
not be applicable (e.g., Chiang et al. 2000). Pos- 
sibly, some of the red wings observed in Seyferts 
are due to a continuum complexity not understood 
at the present time rather than due to an irra- 
diated accretion disk. On the other hand, the 
irradiated disk model, including the associated 
edge/reflection component, works well in some 
other cases, e.g., in the archetypical Seyfert 1 ob- 
ject IC 4329A (Done et al. 2000). In any case, 
a failure of the irradiated disk model cannot be 
considered a justification for fitting X-ray spectra 
with the disk line only and without the associated 
rcfiection component. Instead, new physical mod- 
els for the red wings should be looked for. 

The photoionization process preceding the Fe 
Ka fluorescence should be neglected only when an 
edgc/rcfiection component is so weak that it is ob- 
servationally negligible. This is the case, in partic- 
ular, for lines emitted by a Thomson-thin plasma 
(e.g., Makishima 1986). Such a narrow line is al- 
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ready included in our model, and it is due to both 
extended emission and the emission of the absorb- 
ing clouds. However, such lines are highly unlikely 
to be strongly Dopplcr/gravity broadened as an 
optically-thin accretion flow is necessarily very hot 
close to the central black hole (e.g., Narayan & Yi 
1995; Abramowicz et al. 1995). 

Accordingly, we replace the power-law contin- 
uum in the last model in §3.1 by the sum of a 
power law and Compton reflection normalized in 
such way that the line equivalent width with re- 
spect to the Compton-reflected component alone 
at 6.4 keV is ~ 1.3 keV (George & Fabian 1991). 
This coefficient corresponds to a neutral or moder- 
ately ionized reflection, and we check this assump- 
tion a posteriori below. The reflection component 
(Magdziarz & Zdziarski 1995) is then Gaussian- 
broadened and redshifted in the same way as the 
line. The origin of the Gaussian broadening may 
be the chaotic motion of clouds in the vicinity of 
the black hole. This model yields a much softer 
intrinsic continuum, with F = 1.871^q "q, but very 
strong reflection, fi/27r = SAtfl (Wb ~ 390 cV), 
at x^/^ = 715/637. We assumed here a face-on in- 
clination; for a larger inclination, the fitted f2/27r 
becomes even larger. 

As noted in §3.1, Compton reflection as strong 
as this is unusual for Seyferts. Also, it is likely that 
accretion proceeds via a disk rather than clouds, 
and thus we consider the above model rather im- 
plausible. Therefore we replace the broad Gaus- 
sian line by a line from an accretion disk in the 
Schwarzschild metric with the rest energy of 6.4 
keV and the emissivity oc r~^, extending from an 
outer radius, rout (assumed here to equal lO"') to 
an inner one, rin > 6, where r is in the units of 
GM/c^ (Fabian et al. 1989). The strength of the 
line is coupled to that of Compton reflection in 
the same way as above. The reflection spectrum 
is again relativistically broadened in the same way 
as the line (ref sch in XSPEC). Since we include 
reflection, we constrain the inclination to i > 15°, 
which corresponds to the range at which the reflec- 
tion Green's functions of Magdziarz & Zdziarski 
(1995) are valid. 

The importance of including Compton reflec- 
tion in models of reprocessing features is illus- 
trated in Figure 3 for two cases with strong rel- 
ativistic broadening in an accretion disk. We see 
that the eff'ect is rather strong, and the reflection 



component is by no means negligible. Also, taking 
into account the relativistic broadening of the re- 
flection component is crucial. Roughly, the main 
effect of including relativistically-smeared reflec- 
tion is an overall hardening of the total continuum, 
but without appearance of any distinct Fe K edge. 
The corresponding effect on the fitted parameters 
will be an increase of the fitted F of the incident 
spectrum. 

We also include the Kft component (at the 7.06 
keV rest energy and with the flux of 0.12 of that 
of Ka) for both the disk and the narrow lines. 
Also, we use the optically-thin plasma emission 
with the standard abundances (Mewe et al. 1985) 
instead of bremsstrahlung. This model yields 
X^/u = 701/636, r = 1.63l^;^7, ^f^^ = l-34to:29 
(Wb = 160 eV). The fit strongly prefers an in- 
clination of the disk close to face-on, i = 15+^°, 
and emission of the reprocessed component con- 
centrated close to the minimum stable orbit. Tin = 
6.2tJi, p = 4.0t J;t 

We note here that the quoted equivalent widths 

are given by XSPEC with respect to the contin- 
uum at the peak of the line. E.g., for models 
with /? ~ 4.5, i = 15°, r-m = 6, this peak is at 
~ 5.5 keV (Fig. 3b), at which the continuum pho- 
ton flux is typically by ~ 30% higher than that 
at 6.4 keV. Thus, for the same photon flux of the 
line, the measured equivalent width will be corre- 
spondingly lower. This effect has to be taken into 
account when making comparison with theoretical 
results giving the line strength expected theoret- 
ically, which usually are expressed as the equiva- 
lent width of the line in the local rest frame (e.g., 
George & Fabian 1991; Zycki & Czerny 1994). 

3.3. Complex absorption 

The complexity of absorption in NGC 4151 has 
long been appreciated (e.g. W94). Among others, 
Warwick et al. (1995) and George et al. (1998) 
have c;onsidered absorption by a partially ionized 
medium. However, W99 have found that this 
model provides a bad fit to the present data. 
On the other hand, ionization of the completely- 
covering component in the dual absorber model 
has been found to be required by Piro et al. (2002) 
for BeppoSAX data. We confirm their conclusion 
for our data, obtaining a reduction of by 20 
(the chance probability of 2 x 10~^), a low ion- 
ization parameter of ^ ~ 0.6 erg s~^ cm, but only 
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Fig. 3. — The effect of Compton reflection on the Fe Ka disk hne profile as illustrated by the ratio of the 
total spectrum to the incident power law with T = 1.7. The dashed curve shows the case with reflection 
neglected, the dotted curves show the line profile together with the associated reflection but neglecting its 
own relativistic smearing, and the solid curves give the correct spectrum in which both the line and the 
reflection components are relativistically smeared in the same way. The parameters are Tin — 6, rout — 10^, 
i — 15°, n/2TT — 1 and (a) /? = 3 and (b) /3 = 4.5. The wiggles in the line profile are numerical noise in the 
XSPEC implementation of the disk line in Schwarzschild metric (diskline, Fabian et al. 1989). We note 
that this numerical noise will be averaged out in data from ASCA given its resolution, but it prevents the 
diskline model from using in fits of high resolution data, e.g., those from Chandra. 



small changes in the other fitted parameters. Here, 
^ = AnFion/n, where Fjon is the ionizing flux, n is 
the reflector density, and the temperature of the 
reflecting medium is assumed to be lO'^ K (Done 
et al. 1992). 

Next, we consider possible modifications of the 
partial covering component. Specifically, since 
there is no reason a priori that there should be a 
single partially-covering cloud in the line of sight, 
we investigate the possibility of multiple clouds in 
two ways. First, we replace the partial covering 
component by a Poisson distribution of clouds in 
the line of sight, first proposed by Yaqoob & War- 
wick (1991). However, this yields only a small in- 
crease of r and almost no change in x^; thus, it is 
not required statistically. Furthermore, the Pois- 
son distribution requires an individual event to 
have small probability, and this distribution would 
not be appropriate if there are only a few clouds 
covering a substantial fraction of the source. Thus, 
we consider the case of just two partially covering 
clouds (with the column iV2 covering a fraction /2 
for the second one). This leads to our overall best 
model of the ASCA data, with xV^^ = 655/633. 
The probability that the fit improvement with re- 
spect to the previous model (the dual absorber 



with a weakly ionized full screen, see above) is by 
chance is 4 x 10"''. This model yields i = I7tl°, 
Tin = 6.31q'3, (3 = 5^2' ^ much softer con- 
tinuum, r = l.SltoAo, f^/27r = 1.20t|Jj7 corre- 
sponding to Wt, — 110 eV at the best fit. Inter- 
estingly, these F and fl are close to the Seyfert 
average (Zdziarski et al. 1999). The residuals for 
this model are shown in Figure 2b. 

3.4. The broad-band 0.4-400 keV spec- 
trum 

As noted in §1, the ASCA observation was 
nearly contemporaneous with an observation by 
OSSE, made 1995 Apr. 25-May 9. J97 has found 
that the variability of NGC4151 in soft 7-rays is 
rather slow and modest. Thus, we fit the OSSE 
spectrum (at a free relative normalization) to- 
gether with the ASCA spectra. We find that the 
last model in §3.3 above predicts quite well the am- 
plitude (within 10%) of the OSSE spectrum. The 
model parameters are almost identical as for the 
model of the ASCA data alone (e.g., = 120 cV) 
except that now the e-folding energy is fitted. The 
model P in Table 1 shows the parameters (here, 
the standard values of /3 = 3 and r-^ = 6 have 
been assumed). For this model, we also study the 
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Fig. 4.— The ^S'C^-OSSE spectrum (black, mid- 
dle), compared to the historical range (blue to 
red) of X-ray and soft 7-ray fluxes (highest and 
lowest X-rays from Ginga and EXOSAT, respec- 
tively, Z96; highest and lowest soft 7-rays from 
OSSE, J97, and BeppoS AX /PBS, P.-O. Petrucci, 
private communication, respectively). The solid 
curve gives the model with thermal Comptoniza- 
tion (§3.3). 

allowed ranges of some parameters we otherwise 
keep fixed. We find the relative Fe abundance of 
0.9j^o'3, and, importantly, the reflecting medium 
to be at most weakly ionized, with ^rofl = 1^1° 
erg cm, which range corresponds to the dom- 
inance of Fe i-xii, where effects of reflector's ion- 
ization are negligible. In accordance with this re- 
sult, increasing the disk line energy to 6.7 keV (Fe 
xxv-xxvi) leads to Ax^ = +18, where we made 
conservative assumptions of a free i and inde- 
pendent of Wh- 

Next, we investigate a more physical model 
for the high energy spectrum. The high-energy 
cutoff of the OSSE data is strongly indicative 
of thermal Compton scattering, and the average 
OSSE spectrum of NGC4151 can be well mod- 
elled by it (J97). An e-folded power law is only a 
rough approximation to spectra from that physi- 
cal process. Thus, we now fit the spectrum with 
an isotropic Comptonization model (compps at 
ftp://ftp.astro.su.se/pub/juri/XSPEC/COMPPS) 
of Poutanen & Svensson (1996). The parame- 
ters of the thermal Compton continuum are the 
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Fig. 5. — The components of the model with 
thermal Comptonization (solid curve). The short 
dashes show the intrinsic spectrum of the model 
with thermal Comptonization before absorption. 
The long dashes and dots give the Comptonization 
continuum and the reflection component together 
with the broad Fe K lines, respectively. 

plasma temperature, kT, and the Compton pa- 
rameter, y = ArkT/nieC^, where r is the plasma 
Thomson optical depth. The assumption that 
the seed photons are blackbody at a temperature 
of 5 eV affects the resulting spectrum only very 
weakly. We find the absorber and the soft com- 
ponent parameters are virtually unchanged with 
respect to the previous model, but the reflection 
and disk line (with Wh = 70 eV) are somewhat 
weaker, see Table 1, model C. The normalization 
of the OSSE model with respect to that of SISO is 
close to unity (1.08), and the 2-10 keV index of 
the Comptonization component alone is F = 1.89. 
The spectrum is shown in Figure 4, together with 
the range of X-ray and soft 7-ray spectra histor- 
ically observed from NGC4151. Figure 5 shows 
the main components of the model. 

For this model, we also consider the effect of al- 
lowing for variable (3 and r-m- We find that adding 
these two free parameters results only in a statis- 
tistically insignificant reduction of x^, by —3. At 
this best fit, /3 = 3.9, = 6.3, ~ 60 eV. 

We also note that we find a bad fit of the ther- 
mal Compton model when we use the standard 
dual absorber model for the absorption {x^ /v = 
759/661). Also, that model requires that the nor- 
malization of the OSSE spectrum is up by a fac- 
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tor of 2 with respect to that of the ASCA spectra. 
Thus, inclusion of the OSSE data provides fur- 
ther very strong support for the presence of com- 
plex absorption (with two partial-covering compo- 
nents). This result is in agreement with the OSSE 
spectra from NGC 4151 being consistent with ther- 
mal Comptonization with T ~ 1.9 in X-rays (J97; 
Zdziarski et al. 2000). 

The total intrinsic model luminosity is Lboi = 
3.9 X lO^'^ erg s^^, assuming isotropy and a dis- 
tance of 13 Mpc. The emission measure of the 
optically-thin plasma is then 1.0 x 10^^ cm~^, 
although we stress that we have approximated 
its emission as coUisionally-ionized only (plus the 
0.89-keV line) whereas Ogle et al. (2000) find 
evidence for both coUisional ionization and pho- 
toionization. Our final models have the best-fit 
scattered fraction of fg^ — 1-5% (see Table 1), 
which agrees with theoretical expectations of uni- 
fied AGN models (Antonucci 1993). On the other 
hand, the corresponding model with the standard 
dual absorber has a higher value of fsc — 3.5%. 

3.5. Comparison with W99 

W99 have studied a relatively small subset of 
the data considered here, namely the SIS data 
in the 1-10 keV range. Their model consists of 
the dual absorber, a power law component, Fe 
K disk and narrow lines from a neutral medium 
with both Ka and KP components, and an unab- 
sorbed (scattered) power law. That last compo- 
nent is sufficient to model the soft excess in the 
considered energy range as most of its complexi- 
ties appear below 1 keV (see §3.1), but we note it 
yields a rather high f^c ~ 5%. W99 neglected re- 
flection of the incident power law from the disk as 
they claimed there is no statistical evidence for its 
presence. We note, however, that this makes their 
model only a phenomenological description of the 
data, since their model assumes irradiation of a 
nearly-neutral, optically-thick, disk by the power- 
law source, and under these conditions a part of 
the irradiating photons will be Compton backscat- 
tered (see §3.2). 

We have applied the same model to the SIS 1- 
10 keV data, and have closely reproduced their 
results of r ~ 1.30 and VFb ^ 280 eV at the best 
fit (with x^/^ = 409/382 with our binning and 
data reduction). At this minimum, the width 
of the narrow line component is close to null. How- 
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Fig. 6. — The contour plots (1, 2, and 3a) of the 
photon flux in the disk line, /b, vs. the width of 
the narrow line component, an in the model with 
disk line only (W99), see §3.5. There exist two 
solutions with identical x^j one with cTn ~ and 
a strong disk line, and one with an ^ 0.1 keV and 
the disk line twice as weak. W99 have identified 
only the former. 

ever, we find there exists another minimum with 
the same x^, at which Wb — 150 eV and the width 
of the narrow component is (7^ ~ 0.1 keV, as illus- 
trated in Figure 6. The uncertainties on Wy, given 
by W99, ±20 eV, appear to have been obtained 
under the assumption of an = 0. However, as we 
discuss in §5 below, the presence of strong vari- 
able absorption in NGC 4151 implies the presence 
of Fe K emission of the absorbing clouds (Mak- 
ishima 1986). If the clouds are at a plausible dis- 
tance of ~ lO^GM/c^ and cover > 10% of the 4tt 
solid angle (see §5), Cin will be just '^0.1 keV. 
Thus, the data used by W99 do allow the disk line 
to be substantially weaker than they claim. 

On the other hand, the above ambiguity of 
spectral solutions disappears when the GIS 1-10 
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kcV data arc; added. In that case, there is only 
one minimum, which has Wb — 150 eV and 
(Tn ~ 0.1 keV (similar to other models with disk 
line presented in this work), and setting = 0.01 
keV results in Ax^ = +21. Furthermore, when 
we add the data below 1 keV and model the soft 
excess in the same way as described in §3.1, we ob- 
tain x^/i" = 720/636 (at T = 1.47io;o4, Wb = 110 
eV), which is higher by 19 for the same num- 
ber of degrees of freedom as in the corresponding 
model with reflection and dual absorber described 
in §3.2. Thus, the reprocessed/reflected compo- 
nent associated with the Fe K lines is not only 
needed physically but also it is required statisti- 
cally for our data. 

4. Variability 

Overall, the total flux during the ASCA obser- 
vation varied by ~ 50%. The complexity of the 
spectrum indicated above was born out in the con- 
siderable spectral variability observed. The excess 
variance as a function of energy is shown in Figure 
7. Below 1.5 keV the flux was constant, consistent 
with emission from a extended region. We note 
that, given the energy dependence of the variance 
shown in Figure 7, the average ASCA variance 
of the source is not a suitable quantity for com- 
parison with other, much less absorbed Seyferts 
(Nandra et al. 1997). On the other hand, the ex- 
cess variance in the 7-10 keV range, 7.5 x 10~^, is 
most likely to represent the intrinsic source vari- 
ability (although variability due to changes in the 
absorption cannot be ruled out completely). We 
have found this value to be fully consistent with 
variance of other broad- line Seyferts having intrin- 
sic 2-10 keV luminosity of 1.4 x 10^^ erg s~^ (e.g. 
Leighly 1999), as illustrated in Figure 8. This 
tends to rule out the Compton mirror model pro- 
posed by Poutanen et al. (1996), which would pre- 
dict a reduced amplitude of variability unless the 
mirror has the same size scale as the source. 

Detailed analysis of the spectral variability is 
beyond the scope of this paper. However, we 
note that the highest amplitude of variability is 
in the 2.5-3.0 keV band. This increase with re- 
spect to the 7-10 keV one (likely to be mostly 
intrinsic) provides evidence for the importance of 
variability in the covering fraction/column den- 
sity of the absorber in addition to the spectral in- 
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Fig. 7. — The fractional amplitude of intrinsic 
variability squared as a function of energy. The 
dotted symbol in the 5.5-7.0 keV band accounts 
for the constant narrow component of the emission 
line. 

dex/intrinsic flux variability. We have confirmed 
that our physical model with complex absorption 
and disk line coupled to relativistically-broadened 
reflection also fits the time-resolved spectra. Thus, 
our fit results presented in §3 are not an artifact 
of time averaging. 

The amplitude of variability is reduced in the 
5.5-7.0 keV band compared with neighboring en- 
ergies, as shown in Figure 7. This is most likely 
due to the contribution of the constant narrow 
component of the iron line, which contributes 
about 7% of the flux in that band, implying that 
the excess variance is a factor of 1.15 too low. Ac- 
counting for that makes excess variance in the 5.5- 
7.0 keV band comparable to that in the 7.0-8.0 
keV band. The constant narrow component sup- 
ports our two-component (narrow and disk) model 
of the iron line. The constancy of the narrow com- 
ponent of the line in this data set was first reported 
by Leighly et al. (1998). 

5. The Narrow Line Component 

The flux of the narrow Ka component. In ~ 
2.5 X 10~^ cm~^ s~^ in our best models (see Ta- 
ble 1; somewhat more in our other models), is 
higher than the flux of 1.8 x lO"** cm~^ ob- 
served by Chandra in 2000 March (Ogle et al. 
2000). Ogle et al. find that 65% of the Hne (i.e., 
~ 1.2 x lO"** cm~^ s~^) comes from kpc-scale gas. 
However, a part (dependent on the covering factor 
of the absorber) of the narrow line has to come 
from the absorbing medium (Makishima 1986), 
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Fig. 8. — Comparison of tlie intrinsic X-ray vari- 
ance of NGC4151 with those of other broad and 
narrow- line Seyfert Is (based on Leighly 1999). 
The intrinsic variance of NGC4151 is rather typ- 
ical for other broad-line Seyferts. Note that the 
decrease of the measured variance with increas- 
ing luminosity, L, appears to be an effect of the 
decreasing ratio of the typical length of an obser- 
vation (limiting the range of timescales over which 
the variance is measured) to GM/c^ at M oc L. 
The data shown here are compatible with the total 
variance not showing any trend with L within each 
of the two types of Seyferts (e.g., Leighly 1999). 

which medium is variable on long time scales (e.g., 
Yaqoob et al. 1993). We note that the optical flux 
peaked in 1995-1996, the U band magnitude had 
decreased by AMy w 1.5 by 2000 (Doroshenko et 
al. 2001), and the Chandra observation was made 
when NGC 4151 was in a low state. We suggest 
that the flux of the narrow component was higher 
during 1995 than it was in 2000 (consistent with a 
long-term variability of the soft excess, W94), and 
our results are consistent with those from Chan- 
dra. This is, in fact, confirmed by the findings 
of Weaver et al. (2001), who obtained /„ ^ (3- 
4) X 10"" cm-2 in their fits to five ASCA 
observations of NGC 4151 (in X-ray high states). 
Also, the recent XMM-Newton result of Schurch 



et al. (2002), who found /„ ~ 1.3 x 10"'' cm'^ 
s"^ in a very low state of NGC 4151 (similar to 
the historically lowest state shown in Fig. 4) one 
year after the Chandra observation supports the 
picture of a part of the narrow line responding to 
changes of the intrinsic continuum. 

The narrow component of the Fe K emission 
shows evidence of broadening, with cr^ — 0.05- 
0.1 keV in our best models (see Table 1), and 
somewhat more in disk line models with the stan- 
dard dual absorber in §3.2. This width is re- 
lated to the velocity dispersion at the FWHM, 
Av = (81n2)i/2(cr„/£;K„)c (< 10* km s^^ at the 
above an)- Roughly, Av will be similar to the 
Keplerian velocity, which implies a characteristic 
radius where the emitting clouds are located of 
i?c ~ 2 X lO^GM/c^. These clouds are most likely 
the same clouds that absorb the intrinsic spec- 
trum. As mentioned above, Chandra and XMM- 
Newton measured much narrower Fe Ka lines at 
lower fluxes in later low states of NGC 4151 (Ogle 
et al. 2000). Thus, the "narrow" Fe Ka line ap- 
pears itself to consist of two components: one with 
(Tn ^ 0.1 keV originating in an extended region 
and constant on long time scales, and a compo- 
nent with (Tn 0.1 keV emitted by the absorbing 
clouds and responding to variability of the nuclear 
source. However, given the ASCA resolution, us- 
ing one Gaussian to represent both components of 
the narrow line is sufficient. 

The characteristic line response time to changes 
in the nuclear continuum will be ~ 2Rc/c. For a 
fiducial M = lO'^M©, this gives 2 x 10^^ s, con- 
sistent with the constancy (within the measure- 
ment errors) of the narrow line component during 
our observation (Leighly et al. 1998; see §4 above). 
On the other hand, variability of the absorber dur- 
ing the observation (~ 10^ s) implies, assuming a 
cloud moving in the line of sight with the Keple- 
rian velocity, the characteristic cloud size of a few 
times 10^"^ cm. This size should be similar to the 
size of the intrinsic X-ray source in order to allow it 
to partially cover our line of sight. This, in turn, is 
consistent with the size of the region where most of 
the gravitational energy is dissipated, ~ 20GM/c^ 
or so. 

The equivalent width of the narrow component 
with the parameters of Table 1 is ~ 60 eV with 
respect to other components of the observed ab- 
sorbed spectrum. As calculated by Makishima 
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(1986), the equivalent width from an absorber 
completely surrounding a nuclear source would be 
~ 300 eV at TVh ~ 3 X 10^3 cm-^ (as in Table 1). If 
a half of the observed /„ comes from the absorber 
(and the rest from the extended gas), ~ 10% or 
so of the lines of sight as seen from the nucleus 
should be covered by the clouds. 

6. Discussion and Conclusions 

We have performed a detailed and careful 
analysis of contemporaneous long observations 
of NGC4151 by the ASCA and CGRO satel- 
lites in 1995. We find that proper inclusion of 
the Compton reflection component (taking into 
account its own relativistic smearing) associated 
with the broad Fe K line in the source is important 
for determination of the spectral index, F, of the 
power law component in this source. This com- 
ponent is required physically as well it is required 
statistically for our data (its inclusion results in a 
reduction of by —20 without changing the num- 
ber of free parameters). With it, T > 1.6, which is 
within the 2a range measured for Scyfcrts (Nan- 
dra & Pounds 1994). On the other hand, including 
Compton reflection has a little effect on the fitted 
parameters of the broad line. 

Furthermore, we have studied the complexity of 
the absorber in this source. As previously noted 
by W94, the measured value of T depends crucially 
on the structure of the absorber. So far, the pre- 
ferred model of the absorber has consisted of a uni- 
form neutral screen and a partially-covering cloud 
(the dual absorber, e.g., W94). We have consid- 
ered some simple and plausible modifications to 
this model. First, we confirm a finding by Piro et 
al. (2002) that the uniform screen is weakly ion- 
ized. Furthermore, we find very strong statisti- 
cal evidence (~ 1 — 4 x 10~^) for the presence of 
more than one partial-covering cloud in the line 
of sight. We find that the presence of the second 
absorbing cloud is not only strongly preferred sta- 
tistically, but it also changes the fitted value of F 
to ~ 1.9, close to the average index of Seyferts. 
In this model, the broad line component becomes 
relatively weak, with Wi, - 10^ eV (60-70 cV in 
our preferred physical model with thermal Comp- 
tonization) . The corresponding reflection strength 
is 0/27r ^ 1, also average for Scyfcrts (Zdziarski ct 
al. 1999). Our data also indicate that the broad 



line component is dominated by the Fe K emis- 
sion very close to the central black hole, similarly 
to the case of MCG -6-30-15 (Wilms et al. 2001). 

Obviously, our model is not unique and, e.g., 
further complexity of the absorber is possible. One 
theoretical uncertainty involved regards the equiv- 
alent width of the Fe Ka line with respect to the 
reflection continuum. Our choice of that quantity 
(1.3 keV) is based on George & Fabian (1991), 
whereas Zycki & Czerny (1994) obtain values gen- 
erally lower than 1 keV. Since our fits are driven by 
the line strength, the above choice is conservative 
in the sense that it gives the lowest theoretically 
predicted reflection strength. Adopting results of 
Zycki & Czerny (1994) would yield values of O and 
F somewhat larger than those in Table 1. 

Our model is consistent with a contempora- 
neous observation by the CGRO/OSSE. When 
the joint data are fitted by thermal Comptoniza- 
tion, the plasma temperature is ~ 70 keV, in the 
range of typical values for Seyferts (Zdziarski et 
al. 2000). We have also fitted data from Beppo- 
SAX observations of NGC4151 and have foimd 
excellent fits with the thermal-Compton model 
with two absorbing clouds and Compton reflection 
(work in preparation). 

The flnding that NGC4151 appears to be an 
intrinsically average Seyfert 1 is also conflrmed by 
comparison of its X-ray variance with that of other 
(optically) broad-line Seyferts. The crucial impor- 
tance of complex absorption for the appearance of 
the X-ray spectra is further supported by the en- 
ergy dependence of the variance, which peaks in 
the range where absorption of the intrinsic power 
law is the dominant effect. 

We also have studied the narrow component 
of the Fe K emission. We flnd the presence of a 
broadened Gaussian component with an ^ 0.1 keV 
or so in addition to the broad disk line. This allows 
us to construct a self-consistent scenario (§5) with 
Fe K emission of the absorbing clouds at a char- 
acteristic radius of a few thousands GM/c^. This 
moderately-broad component responds to changes 
of the nuclear continuum and thus was nearly ab- 
sent during Chandra observation of NGC4151 in 
a low state (Ogle et al. 2000). 
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Table 1 

Rest-frame model parameters of the broad-band spectrum 
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Note. — Model P is a phenomenological model with an e-folded power law while model C is a physical Comptonization model. AT, the ionization parameter ^, /n, and {kT, Ec-, 
En, <Tn) are in units of 10^^ cm~^, erg s~^ cm, lO"*^ cm~^ s~^, and keV, respectively. 



